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INTRODUCTION 


It  is  highly  desirable  that  the  motions  of  the  Navy's  modern 
submarines  be  predicted  in  advance  of  full-scale  trials  and  operations 
to  establish  their  safe  operating  envelope  and  their  ability  to  perform 
specific  maneuvers  effectively.  To  predict  these  motions  and  to 
establish  valid  control  strategies  it  is  necessary  both  to  develop  an 
accurate  mathematical  model  of  the  submarine  and  to  determine  accurate 
values  of  the  hydrodynamic  forces  and  moments  acting  on  the  submarine 
hull  and  appendages  which  are  required  in  the  mathematical  model. 

The  David  W.  Taylor  Naval  Ship  R&D  Center  (DTNSRDC)  provided  a 
standard  a?t  of  equations  of  motion  for  use  in  submarine  computer  simula- 
tion in  Reference  1.  These  equations  have  been  used  to  simulate  the 
trajectories  and  responses  of  submarines  in  six  degrees  of  freedom 
resulting  from  various  types  of  normal  maneuvers  as  well  as  for  extreme 
maneuvers  such  as  those  associated  with  emergency  recoveries  from  a 
stamp  lane  Jam. 

Reference  2 gives  a general  description  of  the  effort  at  the  Center 
to  predict,  evaluate,  and  improve  the  stability,  control,  and  maneuver- 
ing characteristics  of  the  Navy's  submarines.  Including  modifications 
and  improvements  made  to  the  equations  of  motions  given  in  Reference  1. 
The  Improvements  to  the  equations  of  motion  outlined  in  Reference  2 have 
resulted  in  better  correlation  with  full-scale  trial  data. 

This  report  has  been  prepared  to  provide  those  working  in  the  field 
of  submarine  stability,  control,  and  maneuvering  with  documentation  of 
the  current  interim  mathematical  model.  This  report  defines  the 
notation,  axos  systems,  and  sign  conventions  used  for  the  equations 
and  presents  the  DTNSRDC  Revised  Standard  Submarine  Equations  of  Motion 
for  performing  computer  simulation. 


NOTATION 


*i‘  *V  ci 

£ 

- 

Sets  of  constants  used  in  the  repre- 
sentation of  combined  thrust  and 
drag*  in  the  axial  equation 

A 

■■  y 

v-$ 

Projected  area  of  hull  plus  deck 
in  xx-plane 

A 

X 

A 

A ' | 

Z i2 

Projected  area  of  hull  in  xy-plane 

b(x) 

to’  - 

Local  beam  of  hull  in  xy-plane 

J b(y)dx  - Az 

it 

B 

B*  " r ---2-r 

Y pt  V 

Buoyancy  force  of  envelope  displace- 
ment , positive  upward 

C 

Variable  coefficient  used  in  scaling 
model  thrust  and  drag  data  to  full- 
rcale.  Function  of  AX 

V C7*  C8 

Constants  used  in  computing  C 

CB 

Center  of  buoyancy  of  submarine 

Cd 

c - CFD 

d i^v2 

Coefficient  used  in  integrating 
forces  and  moments  along  hull 
due  to  local  cross-flow 

CFD 

Cross-flow  drag 

CG 

Center  of  mass  of  submarine 

Modified  nondimensional  sectional 
lift-curve  slope  used  in  computing 
the  effects  of  the  hull-bound 
vortex  due  to  lift  on  the  bridge 
fairwater 

h(x) 

h(x)  * - 

Local  height  of  hull  plus  deck  in 

xz-plane  j h(x)dx  - Ay 

£ 


xy 


y* 


IX 


Kj 


X 


PlP 


Homdc  of  Inertia  of  submarine 
•hout  x axis 


Movent  of  inertia  of  submarine 
about  y axta 


Movent  of  inertia  of  submarine 
about  i axle 


Product  of  inertia  with  respect 
to  the  x and  y axca 


Product  of  inertia  with  respect 
to  the  y and  x axes 


Product  of  inertia  with  reapect 
to  the  s and  x axes 

Hydrodynamic  movent  component  about 
x axis  (rolllug  moment) 


Coefficient  used  in  representing 

2 

K as  a function  of  u 


Coefficient  used  in  representing 
K due  to  interference  effects  of 
vortices  from  the  bridge  fairwater 
on  the  stern  coutrol  surfaces 


Coefficient  used  In  representing 
K as  a function  of  up.  Does  not 
include  effects  of  vortices  from 
the  bridge  fairwater  on  the  stern 
control  surfaces 


Coefficient  used  in  representing 
K as  a function  of  £ 


6 


Coefficient  used  In  representing 
K as  a function  of  p|p| 


Coefficient  used  tn  representing 
it  as  e function  of  qr 


Coefficient  used  In  representing 
K as  a function  of  ur.  Does  not 
Include  effects  of  vortices  fro® 
bridge  fairvater  on  stern  control 

surfaces 


Coefficient  used  In  representing 
It  it  i function  of  t 


Coefficient  used  In  representing 
K a*  s function  of  uv.  Does  not 
include  effects  of  vortices  from 
bridge  fairvater  on  stern  control 
surfaces 


Coefficient  used  In  representing 
K «s  s function  of  v 


Coefficient  used  In  representing 
K as  s function  of  wp 


Coefficient  used  in  representing 

2 , 

It  in  i function  of  u 6 

r 


Coefficient  used  in  representing 
K as  a function  of  (n  - ^;C 

Coefficient  used  In  representing 
K due  to  at  the  stern  control 
surfaces 


Coefficient  used  In  representing 
K due  to  et  the  stern  control 
surfaces  ' 

Overall  length  of  submarine 


JO*3 


5 piV 


■ r-T 

5 


M ' • 

q 1 .4 

2 p' 


M, ' • a~ 

A 1 .5 

5 pi 


M • - 

fp  i Pi' 


H ’ - r~V 
W \ Pi3 


M • * 

* 3 P«A 

"M'-fS- 


H . . ■ . -gJyl* 
"|w|*  l p«3 


M ' . W 

W \ Pi3 


M • . J™. 

6b  ipe3 


M . - 

63  i Pi3 


Hats  of  submarine , Including  water 
in  f ree-f looding  spaces 

% 

Hydrodynamic  moment  component 
about  y axis  (pitching  moment) 

Coefficient  used  in  representing  H 

2 

as  a function  of  u 

Coefficient  used  in  representing 
H as  a function  of  uq 


Coefficient  used  in  representing 
M as  a function  of  $ 


Coefficient  used  in  representing 
M as  a function  of  rp 


Coefficient  used  in  representing 
M as  a function  of  uw 


Coefficient  used  in  representing 
M as  a function  of  w 


Coefficient  used  in  representing 
M as  a function  of  u|w| 

Coefficient  used  in  representing 
M as  a function  of  w|(v^  ♦ 

Coefficient  used  in  representing 
M as  a function  of  |w(v^  + 

Coefficient  used  in  representing 

2 

M as  a function  of  u 6. 

D 

Coefficient  used  in  representing 

2 

M as  a function  of  u 6 

s 


Hydrodynamic  moment  component  about 
t axle  (yawing  moment) 


Coefficient  uaed  In  representing 

2 

N aa  a function  of  u 


Coefficient  used  In  representing 
N aa  a function  of  up 


Coefficient  used  in  representing 
N as  a function  of  ^ 

Coefficient  used  in  representing 
N as  a function  of  pq 


Coefficient  uaed  in  representing 
N mr  a function  of  ur 


Coefficient  used  In  representing 
N as  a function  of  r 


Coefficient  used  in  representing 
N as  a function  of  uv 


Coefficient  used  in  representing 
N as  a function  of  v 


Coefficient  used  in  representing 
N as  a function  of  v | (v^+  *| 

Coefficient  used  in  representing 

2 

N as  a function  of  u 6 

r 

Coefficient  used  in  representing 
N as  a function  of  ~ C 

Angular  velocity  component  about 
x-axis  relative  to  fluid  (roll) 


AnguVr  acceleration  component 
about  x-axis  relative  to  fluid 

Angular  velocity  component  about 
y-axis  relative  to  fluid  (pitch) 


Angular  acceleration  coaponen 
about  y-axls  relative  to  fluid! 

Contribution  of  propeller  torque 
to  K and  Machinery  equation 


Angular  velocity  coaponent  about 
x-axis  relative  to  fluid  (yaw) 


Angular  acceleration  component 
about  x-axis  relative  to  fluid 

Constanta  used  in  coaputlng  *2 


Tine 

Velocity  of  origin  of  body  axes 
relative  to  fluid 


U ' 

a 


(u2+v 


2.  2 

a 


1/2 


U 


Velocity  of  aternplane  x-coordinate 
relative  to  the  fluid 


Coaponent  of  U in  direction  of  the 
x-axia 


Time  rate  of  change  of  u in 
direction  of  the  x-axis 


u 

£ 

uc*  * — - Comound  apoed:  ateady  value  of 

ahead  speed  component  u for  a 
given  propeller  rpm  when  body 
angles  (a,$)  and  control  surface 
angles  are  xero.  Sign  hanges 
with  propeller  reversal 
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V 


v' 


V 

U 


V 


V 

s 


'FW 


VFW(t  " tW) 


Component  of  U in  direction  of 
the  y-axia 

Time  rate  of  change  of  v in 
direction  of  the  y-axis 


Velocity  component  in  the  y-axis 
direction  at  the  quarter  chord 
of  the  sternplanes.  vg  - v + xgr 


Velocity  component  in  the  y-axis 
direction  at  the  starting  position 
of  the  hull-bound  vortex  due  to 
lift  on  the  bridge  fairwater. 


FW 


+ x^r  - 


ZFWP 


Velocity  component  in  the  y-axis 
direction  at  the  bridge  fairwater. 

VFW  " V + XFWr  ~ ZFWP 

Value  of  at  time  » t - x(x) 


VFW(t  - TX) 


Value  of  at  time  ■ t - 


v(x) 


V(X)' 


Velocity  component  in  the  y-axis 
direction  of  any  x coordinate 
v (x)  - V + xr 


w 


w 


w 

8 


v(x) 


w 


w 

u 


w 


U 


w(x)  ' 


U 


Component  of  U in  the  direction 
of  the  z-axis 

Time  rate  of  change  of  w in 
the  direction  of  the  z-axis 


Velocity  component  in  the  z-axis 
direction  at  the  quarter  chord 
of  the  sternplanes  w * w - x q 

8 S 


11 


Velocity  component  in  the  z-axis 
direction  of  any  x coordinate 
w(x)  ■ w - xq 


Mwm  *sgg*re  ^yrwaqgwywyy  rc^r^ti1 


m*mnm\  mm 


\ ! 

V 

w<  -t-JL  3 

5 o’  u 

Weight  of  aubmarine,  including 
water  in  free  flooding  epacee 

t ; 

► 

i . 

i 

!■ 

| 

X 

, X 

* - J 

Longitudinal  body  axie;  alao 

the  coordinate  of  a point  relative 

to  the  origin  of  the  body  axia 

5 

r 

X1 

* . *L 
xi  * I 

The  x coordinate  of  the  atarting 
noeition  of  the  hull-bound  vortex 
due  to  lift  on  bridge  fairwater 

i 

[•: 

X2 

, x2 
xa  * r 

The  x coordinate  of  the  aft-moat 
noeition  of  the  hull-bound  vortex 
due  to  lift  on  bridge  fairwater 

' 

f 

XB 

, X» 

*B  “ 1 

The  x coordinate  of  the  LB 

6 

i: 

r 

XG 

, XG 
*G  '*  T 

The  x coordinate  of  the  CC. 

r 

X 

o 

X 

, ° 

Xo  *4 

A coordinate  of  the  displacement 
of  the  origin  of  the  body  nxia 
relative  to  the  origin  of  a set 
of  fixed  axes 

* 

XR 

, XR 

V *r 

The  x coordinate  of  the  quarter 
chord  of  the  rudders 

i 

xs 

, XS 

xs  r 

The  x coordinate  of  the  quarter 
chord  of  the  aternplanea 

f 

£'■ 

XT 

xt  " T 

The  x coordinate  of  the  average 
location  of  the  aternplanea  and 

ruddera.  xt  * 5 ^XS  + XR^ 

r 

XAP 

, XAP 

XAP  " l 

The  x coordinate  of  the  after 
perpendicular 

FW 


FW 


FW  * 


The  x coordinate  of  th®  quarter 
chord  of  the  bridge  fairwater 
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i t iftti  jfcilii 


' ■™':F-  Jf^rp  fr^i- - 


X' 


\ piv 


ax 


AXr  AX2,  AX3 
X X ' 

qq  qq 


. Xm - 

X 


Hydrodynamic  force  component, 
along  x-axis  (longitudinal,  or 
axial,  force) 

Variable  coefficient  uaed  in 
scaling  model  thrust  and  drag 
data  to  full  scale 

Constanta  uaed  in  computing  AX 


Coefficient  uaed  in  representing 
X as  a function  of  q*\ 


X 


Coefficient  used  in  representing 
X as  a function  of  rp 


Coefficient  used  in  representing 
X as  a function  of  r 


Xu 


X. 

X, ' - ~ 

U \ P*3 


X 


vr 


X ' 


vr 


vr  i pi3 


vv 


vv 


vv  j pp 


Coefficient  used  in  representing 
X as  a function  of  u 


Coefficient  used  in  representing 
X as  a function  of  vr 


Coefficient  used  in  representing 

2 

X as  a function  of  v 


ww 


X ’ 


ww 


X 

__W9_ 

-X 


Coefficient  uaed  in  representing 
X as  a function  of  wq 


Coefficient  used  in  representing 

2 

X as  a function  of  w 


X6b6b 


y i 

<Sb6b 


<5b<Sb 
| P*2 


Coefficient  uaed  in  representing 

2 2 

X as  a function  of  u <5, 


n 


6r6r 


t < - Mr 

Mr  \ P*2 


os6s 


t * , X|Ss(:)8 

‘«s6»  1 ft#2 

2 p* 


y’-f 


, yB 

V * r 


V • . 

yG  £ 


v « . 
yo  l 


Y’ 


5 p*  V 


V * m i 

* J P*2 


Y 1 «- 

- i p*3 


Y.  • ■ 

f Jp** 

Y 

Y j | ’ ■ 

PIPI  i p*A 


JELSL 


Y ’ «-  . . 

Pq  \ ^ 


,«****«■*' 


lW  IIIJI.WI.I1  " I ■ ■inww 


Coefficient  used  in  representing 

2 2 

X as  a function  of  u 5 


Coefficient  used  in  representing 

2 2 

X as  a function  of  u 6 

8 

Lateral  body  axis;  also  the 
coordinate  of  a point  relative 
to  the  origin  of  body  axes 


The  y coordinate  of  CB 


The  y coordinate  of  CG 


{ r 


A coordinate  of  the  displacement 
of  the  origin  of  the  body  axis 
relative  to  the  origin  of  a set 
of  fixed  axes 

Hydrodynamic  force  component  along 
y axis  (lateral  force) 

Coefficient  used  in  representing 

2 

Y as  a function  of  u 


Coefficient  used  in  representing 
Y as  a function  of  up 


Coefficient  used  in  representing 
Y as  a function  of  £ 


Coefficient  used  In  representing 
Y as  a function  of  pjp] 


Coefficient  used  in  representing 
Y as  a function  of  pq 
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Y 


r 


1 

l 


pi 


3 


R 


v v R 


p£ 


Coefficient  used  in  representing 
Y «b  a function  cf  ur 


Coefficient  used  in  representing 
Y as  « function  of  r 


Coefficient  used  in  representing 
Y as  a function  of  uv 


Coefficient  used  in  representing 
Y as  a function  of  v 


Coefficient  used  in  representing 

. , , i , 2 . 1/2 

Y as  a function  of  v|(v  + w ) 


Coefficient  used  in  representing 
Y as  a function  of  vp 


Coefficient  used  in  representing 

2 

Y as  a function  of  u 6 


Coefficient  used  in  representing 
Y as  a function  of  (n  ~ 


Normal  body  axis;  also  the 
coordinate  of  a point  relative 
to  the  origin  of  the  body  axis 


The  z coordinate  of  the  hull 
centerline 


The  z coordinate  of  the  CB 


The  z coordinate  of  the  CG 


A coordinate  of  the  displacement 
of  the  origin  of  the  body  axle 
relative  to  the  origin  of  a act 
of  fixed  axee 


The  s coordinate  of  the  42>percent 
apan  of  the  bridge  falrvater 

Hydrodynamic  force  component 
along  s-axia  (normal  force) 


Coefficient  uaed  in  representing 

2 

Z as  a function  of  u 


Coefficient  used  in  representing 
T as  a function  of  uq 


Coefficient  used  in  representing 
Z ss  a function  of  4 


Coefficient  used  in  representing 
Z as  a function  of  vp 


Coefficient  used  in  representing 
Z as  a function  of  uw 


Coefficient  used  in  representing 
Z as  a function  of  w 


Coefficient  used  in  representing 
Z as  a function  of  u|t'| 

Coefficient  used  in  representing 
Z as  a function  of  |w(v2  + v2)1^2! 


Coefficient  used  in  representing 

2 

Z as  a function  of  u 6. 

D 


a* 


x (x)  Tin#  interval  required  for  vortex 

to  travel  f roe  x^  to  any  x 

coordinate  eft  of  x^.  It  ie 

implicitly  defined  ae: 
t 

J u(t)dt  - xx  - x 

t-T(x) 

p Maee  density  of  water 
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ayaMfciaiiiiiiMi 


axial  force  equation 


fu  - vr  + VS  - *g(q*  + r*)  + Vgtpej  - t>  + v-c(pr  + 4) J 


♦ ! l*  [V ,J  + V rJ  + S’  tpj 

♦$‘3[V4*V  vr  + V'"'] 

♦ 5 *2  [v  *2 " V “2] 


. .» 


+ | t2  v 


6t6r 

-(W-B)  Bln  0 ♦ F 


4r2  + W “2{.2  + W “V] 


*P 


*P 


T - DRAG 
P 


( £ t2  ^ 4-  AX)  u2  ♦ bt  Cuuc  + ct  C2uc2] 


AX, 


where  AX  ■ AX^  + 


(AX3  + 1o81q  u) 


C - C6  + (C?  + Cg  AX) 


1/2 


n 


Mote: 


F Is  represented  by  T - DRAG  when  propulsion  characteristics 
xp  P 

•re  available.  Otherwise  the  second  expression  for  F Is  used. 
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LATERAL  FORCE  EQUATION 


- up  ♦ ur  - ys(r2  ♦ P2)  ♦ «G(qr  ” ♦ *G(qp  ♦ *)J 

2 |4[V  V pIpI  * V "] 

f t3  [V  ur  ♦ Yp‘  up  ♦ y ♦ ♦ up] 


♦ £ 
♦ fi 


♦ ft2  [y/  u2  ♦ Yv’  uv  4-  Yv|vjR'  v | (v2  4-  U2)l/?|] 


+ i'2  [V  *V  W UMn-0C] 

- f Cd  J h(x)  v(x)  Jlw(x))2  4-  lv(x)]2}  1 


n 


dx 


X1 
» u 


w(x)  VpyCt  - i (x))  dx 


4-  (W  - B)  cost)  aln<F 
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NORMAL  FORCE  EQUATION 


[*  - uq  + vp  - «c(pJ  + + ^(rp  - 4)  ♦ ygtri]  + 4)J  • 

♦ ! **  V « 

+ f *J  [V  * * V » + V *»] 

♦ ft2  [z.’  u2  + Zw' 

1/2  1 

♦ $ *2  [«W-  «W  ♦ Zw'  l"<v2  ♦ w2)  |J 

♦ § *j  [V  * V ♦ h.n  “\  ("  - e)  c] 

-fc,  J b<x>  «<x>  J[»(x)l2  + [v(x)]2J  1/2  dx 

*1 

♦ £ l£L  J v<*>  ^FW(t  • TW>  dx 

*2 

+ (W  - B)  co«9  coa$ 


„* 


"V 


JWWWli 


ROLLING  MOMENT  EQUATION 

lx  1>  + <*,  - Iy)‘»r  ~ (*  + M)I1X  ♦ (r2  “ <|2>Iyg  ♦ <P*  " $>1^ 

♦ ■ [rG(*  “ uq  ♦ vp)  - «G(v  - wp  ♦ ur)j  • 

♦ ft5  [y  t ♦ V t ♦ V v + ^1*  p|p|] 

♦ f »*  [%•  «P  * V « ♦ V * + V «p] 

+ 2 *3  [*»’  "2  * V “»  * V uvPW(t  * tT>] 

+ f ,3[V  “V*W  "2{r  (x-i)*] 

+ f t3<“2  + V + “s2)  Ss2  [V  ‘*s  + V ,ln  8*s] 
*! 

+ 2 t2,l*  CL  j W(x)  vFV(t  “ Tlxl)  dx 
X2 

+ (YqW  - y^B)  co»0  cos$  - (*GW  - *bB)  co«6  »in$ 


I -'■J*-1-  "■  in  1 , ■ jp4»i.u»iwwuwi  m i. . . >, ., w.  u 


PITCHING  MOMENT  EQUATION 


Iy  4 * dx  - It)rp  - (t  + qt)!^  + <P2  - '2)1I«  + <«»  - m 
+ m [,g(4  - vr  + «q)  - *q(»  - »q  + »P)]  ■ 

+ § t5  [V  4 + Mtp-  rp] 

+ ! t‘  [m4'  * + m,'  “q] 

1/2  "I 

+ !j  l3  k'  u2  + Mv’  U»  + MV|V|R’  »|  (V2  + v2)  I J 

+ 2 *3  [M|w| ' u|u|  + ,W  ^ + "2)  1 J 

+ \ *3  [»{,'  “k  + V *\  + M«.n’  "2{«(n  ' e)  c] 


+ 2 c . f X b(x)  w(x)  J[w(x)]2  +[v(x)]2 
2 d 


J. 


1/2 


dx 


- e ACL  J X v(x)  v^t  - T[x])  dx 
*2 


- (xG  H - Xj  B)  cos©  COS*  - (*G  w - Zg  B)  sin© 


YAWING  MOMENT  EQUATION 


1z  * + (Iy  ’ Ix)pq  “ + rp>Iy8  + <q2  - P2)Ixy  + (rq 

+ ■ [*G<v  - wp  + ur)  - yG(u  - vr  + vq)]  - 

+ \ *5  [v  * + y * * v pi] 

+ ^ ^ j^Np'  up  + Nr*  ur  + N^' 

+ \ ^ [V  u2  + V uv  + Nv|v|r'  vl(y2  + w2>1/2| 

♦^[v  «V>W  »2«r(n-0C] 

“ 2 Cd  J x h(x)  v<x>  +[v(x) ]2  | lU  dx 

p . A 

~ 2 ACL  J x w(x)  vFW(t  ' T{xJ)  dx 
x2 

-*■  (XG  W - xB  B)  cos©  Bind)  + (yQ  W - yfi  B)  sin© 
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AUXILIARY  EQUATIONS 


k 


i 

t 


■ • 

♦ ■ p + ain0 

e - q co«<t>  - r sin<J> 

* 

♦ ■ (t  coa$  + q aln$)/coa9 

*o  " u co*®  ‘'‘’•’I'  + v <®in$  *in0  coal l>  - co*p  etni |») 
+ w (aln$  ainif/  + coa$  aind  coaiJO 

9Q  ■ u coa0  alntj/  -f  v (coa$  coal p + ain$  a in 6 ainifO 
+ V (co*4>  aln6  alniji  - ain$  coatfi) 

*o  " “ u ®ln0  + v CO80  + v coa0  co»<p 

U - (u2  + v2  + v2)1/2 

V lor  lei  i 6st 

*2  *' 

*1  ~ <*i  - V <S1  + SjlBl)  tor  | S|  » B 


1 


, ipvp«=r,w^yjr^'fT\r*F*'»s 
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